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Extended Abstract 

A methododlogy for very fast design of 3DOF entry trajectories subject to all common 
inequality and equality constraints is developed. The approach make novel use of the well- 
known quasi-equilibrium glide phenomenon in lifting entry as a center piece for convenirently 
enforcing the inequality constraints which are otherwise difficulty to handle. The algorithm 
is able to generate a complete feasible 3DOF entry trajectory, given the enttry conditions, 
values of constraint parameters, and final conditions in about 2 seconds on a PC. Numerical 
simulations with the X-33 vehicle model for various entry missions to land at Kennedy Space 
Center will be presented. 

1. Entry Guidance Problem Formulation 

The 3DOF equations of motion of the RLV over a spherical Earth in an Earth-centered 
and Earth-fixed coordinate system are given by 1 
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where r the radial distance from the center of the earth to the vehicle. The longitude and 
latitude are 6 and 0, respectively. The Earth-relative velocity is U, and D = pV 2 S re fC D /2m 
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and L = pV 2 S ref C L /2m are aerodynamic drag and lift accelerations, where S Tef is the ref- 
erence area of the RLV, m the mass, and p the atmospheric density. The lift and drag 
coefficients C L and C D are modeled as functions of angle of attack a and Mach number. 
The atmospheric density p can be approximated by an exponential function of the altitude 
p = p oe -/»(r-Ho) j w ith Rq = 6378135 m, 0 > 0 and p 0 > 0 being constants. The Earth-relative 
flight path angle is 7 and o the bank angle. The velocity azimuth angle ip is measured from 
the North in a clockwise direction. Note that while the rotation of the Earth is important 
in simulating the trajectory of the RLV, it may be safely ignored in the entry trajectory 
planning and entry guidance law development. The feedback nature of the closed-loop entry 
guidance will easily compensate for the difference. 

The entry gui dan ce problem is defined as follows: given the entry conditions at the entry 
interface, the coordinates of the landing site, find the required bank angle profile a (t) and 
angle of attack history a(t) such that 

1. The RLV reaches the terminal area energy management (TAEM) interface with speci- 
fied conditions on altitude, velocity, velocity azimuth, and range to the heading align- 
ment cone (HAC); 

2. The trajectory observes path constraints that specify, respectively, maximum heating, 
load/acceleration, dynamic pressure, and equilibrium glide condition; 

3. The a-profile is consistent with the RLV trim capability, and both a and a do not 
exceed the flight control system authority. The latter constrains the maximum magni- 
tude, rate and acceleration of a and a. 

The entry guidance problem is traditionally addressed in two parts: 

(1) off-line design of a reference: This is where some type of reference trajectory/profile is 
iteratively designed to meet the above requirements; 

(2) on-board feedback tracking of the reference and minor adjustments when necessary. 

The capability of generating a feasible reference trajectory on-board will provide great 
flexibility to the RLV operation and eliminate a most labor-intensive part from the entry 
guidance planning. The two key ingredients for a fully automated entry guidance system are 
on-board generation of a reference trajectory and a tracking control law that automatically 
adjust to the trajectory generated (no gain-scheduling needed). The second task has been 
successfully addressed by a recent development.^ The first task is addressed here. 

2. Rapid Generation of 3DOF Entry Trajectories 

The chief challenges of on-board entry trajectory generation lie in designing the 3DOF 
trajectory to meet all the path (inequality) constraints reliably in all cases. Even in off-line 
trajectory planning, the commonly employed methods for trajectory design will inevitably 
require many iterations. Often times it takes experience, good understanding of the method 
used, and insight into the mechanics of the entry flight to ” tweak” the process to get sat- 
isfactory results. In other attempts to generate the entry trajectory on-board (e.g., ), the 
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path constraints are not considered at all. Without the ability to ensure path constraints 
not to be violated, the usefulness of such an algorithm would be very limited. Since the 
dynamics of the RLV in Eqs. (1-6) are highly nonlinear and lack any special features such as 
differential flatness, the success of an on-board trajectory generation algorithm relies on clev- 
erly decomposing this infinite-dimensional search problem into a finite-dimensional search 
problem with a very small dimension. 

2.1 Assumptions and Objective 

We shall refer to a solution of the differential equations (1-6) (state and inputs cr(t) and 
a(t)) subject to all the imposed conditions listed in Section 1 for a given RLV as a feasible 
trajectory for the RLV. The following assumptions establish the scope of applicability of the 
algorithms to be presented. 

1. The entry flight of the RLV is lifting entry ( L/D ^ 0). 

2. The TAEM conditions are specified in terms of altitude, velocity, range-to-HAC, and 
velocity heading pointing to the HAC. 

3. All the path constraints can be expressed as constraints in the altitude and velocity 
space (when a is given). 

4. There exists at least one feasible trajectory for the given entry conditions, imposed 
TAEM conditions, path constraints, and the given vehicle capability. 

5. A nominal a-versus-Mach profile is available, and limited variations about this nominal 
profile are allowable. 

The last assumption is not a necessity for the algorithms, but a practical condition. An 
arbitrarily designed a(t) profile may not be compatible with the requirements due to thermal 
protection and flight trim conditions. And a nominal a profile, once determined for a given 
vehicle, typically does not change significantly from mission to mission. The Assumption 3 
encompasses all the typical constraints considered in entry flight. 

The objective is to find a feasible 3DOF trajectory within a short time period that is 
comparable to on-board guidance environment. Such an algorithm not only can be applied 
on-board, but will also provides a fast tool for off-line trajectory generation and analysis. 
It should be noted that no optimization of any particular performance index is explicitly 
attempted in the design. But for entry flight, the trajectory is highly constrained by the 
path cosntraints and TAEM conditions so that an optimal and a feasible trajectory do not 
differ appreciably. 

2.2 Quasi-Equilibrium Glide Condition 

The centerpiece in this algorithm for path-constraint enforcement is the use of the quasi- 
equilibrium glide condition. It is well known in flight mechanics that along a major portion 
of a lifting entry trajectory, the flight path angle 7 is small and varies relatively slowly. For 
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a constant bank angle cr, setting cos 7 = 1 and 7 = 0 in Eq. (5) produces the so-called 
equilibrium glide condition 

L cos <7 + (V 2 — — )- = 0 (7) 

r r 

In actual flight, a is not a constant, but the above observations are still true. In other words, 
for a major portion of the lifting entry flight, at any point of the trajectory where the pair 
(r(t), V(t)) is given, the value of the bank angle at that point will be close to the correspond- 
ing value of a(t) from condition (7). Clearly the value of a from Eq. (7) is generally different 
from one point to another along the trajectory. Henceforth we shall call the condition (7) the 
quasi-equilibrium glide condition (QEGC) to signify that we allow time-varying a in Eq. (7). 

The QEGC provides us a simple and effective means to shape the altitude and velocity 
profiles. For any given value of V at a point of the trajectory, if a value of <j is selected, 
the value of altitude at that point is determined from the QEGC. If a is selected between 
a maximum and minimum value, the corresponding altitude will be between a minimum 
and maximum range. The obvious lower bound for o in this process is zero degree, or a 
fixed value when the equilibrium glide condition is desired to be enforced at that value. The 
maximum vahie of a is determined in the following way: Suppose that the path constraints 
are 


Si(r, V y a) < Ci, i = l,...,l ( 8 ) 

where q > 0 is constant. The region in the r-V space where all the above conditions are 
met is called the entry flight corridor. For all the commonly used entry path constraints, the 
value of Si(r, V, a) increase as the altitude decreases (atmospheric density increases). Figure 
1 shows the nominal entry trajectory for the X-38 Crew Return Vehicle and several typical 
path constraints. 
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Fig. 1: Entry trajectory for the X-38 Crew Return Vehicle 

On the constraint boundary Si(r, V, a) = c*, for any given V, the value of r is determined 
because a is from the nominal a-versus-Mach profile. This r represents the lowest altitude 
at this velocity without violating the i-th constraint. Find the corresponding bank angle 
from the QEGC, and denote it by < 7 l max (V) because it will be the maximum allowable bank 
angle at this velocity without violating the constraint Si(r , V, a) < Q. Define for all V in the 
portion of the entry trajectory where the QEGC is valid 

<Tmax{V) = min W'maxiY)} (9) 

l<i<t 

Then, for any such V, when we choose a a < <J max {V), the corresponding r from the QEGC 
will ensure that all the path constraints in Eq. (8) are satisfied. Figure 2 shows the variation 
of such (Tmax ( V) for the X-38 and the corresponding boundaries of the entry flight corridor. 
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Fig. 2: X-38 Entry flight corridor and a from QEGC for the path constraints 


2.3 Altitude- and Velocity-versus-Range Profiles 

The QEGC will be used to facilitate the search for r- and F-versus-range profiles sat- 
isfying path and TAEM constraints into a one-parameter search. All the steps described 
hereafter are automated. Note that the QEGC can only be satisfied by some triplet (r, V, <r) 
once sufficient dynamic pressure has been built up. For orbital entry, there will be an initial 
period of the flight where the atmospheric density is too low for the QEGC to be met by 
any triplet (r, V, cr). In such a case in the trajectory design process we start with an initial 
descent phase where a = a o, a constant. 


Initial Descent 

The constant a 0 is chosen such that the initial descent will enter the entry flight corridor 
and transition smoothly onto the QEGC. Since the heading of the RLV at entry interface 
is pointed toward the landing site, the differential equation for the range-to-go to the HAC 
may be written as 


Stogo 


V cosy 


( 10 ) 
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Using Stogo 8 s the independent variable, and for the chosen <7o and the nominal a-profile, the 
Eqs. (1), (4) and (5) are numerically integrated. The QEGC is continuously checked during 
the integration. Once the QEGC begins to admit a solution a for the pair (r, V) obtained, 
we start the checking for transition point onto the QEGC. The most effective criterion for 
determining the transition point we found is to compare the slope of dV/dr on the numeri- 
cally integrated descent portion with that along the QEGC. Whenever the two slopes match 
within a tolerance, the trajectory transitions onto the QEGC. 

The value of <r 0 is first set to be zero. But in some cases <r 0 = 0 will result in the trajec- 
tory bouncing a few times before it can transitions onto the QEGC (characterized by large 
oscillations in the altitude history, a known phenomenon in entry flight). If this undesirable 
bouncing happens, a new oo > 0 is selected in a pre-selected increment and the numerical 
integration restarted from the entry interface. Note that the path constraints are all inactive 
in this phase because the atmospheric density is still very low. At most 3-4 repetitions, 
this process will have been completed, and the trajectory is on the QEGC. The integrated 
r(stog 0 ) and V(sto 9 o) are stored for later use. 

Determination of Bank Angle Profile on QEGC 

The values of a at the first point and the last (the TAEM point) on the QEGC are solved 
Eq. (7), since the values of r and V are already known or specified. Denote them by cri and 
(t T aem- Let a m i d be a value to be determined at the midway point of the TAEM velocity 
and Vi, the velocity at the transition point from the initial descent to QEGC. A piecewise 
linear, continuous a(V) profile can be obtained using <Ti, a mid and (Ttaem (or a quadratic 
<r(V) if an additional condition d<r{V)/dV = 0 at o' mid is used). Again in Eq. (4), use St ogo 
as the independent variable, sin 7 « 0, cos 7 ~ 1, and replace D = {Cd/Cl)L where L is 
substituted from the QEGC (7) 

dV _ /d _ y2\ {CdIC£) ^ 

ds togo \r ) V cosa 

For any selected o m id, <x(V) is determined for all V € \Vtaemi Vi]. At any V and the corre- 
sponding o’ ( V ) , the value of r is found from the QEGC. From the nominal o: profile, and the 
given RLV aerodynamic model, Cd /Cl is evaluated. Thus Eq. (11) can be numerically in- 
tegrated from the s togo at the transition point onto the QEGC to the required range-to-HAC 
at the TAEM interface. The resulting TAEM velocity is obtained. If the resulting TAEM 
velocity is not equal to the required Vtaem > &mid is adjusted and the process repeats. The 
corresponding r profile is also obtained by evaluating the value of r from the QEGC at any 
integrated velocity value. Note that the r-versus-range profile obtained this way satisfies 
the required TAEM altitude condition automatically if the TAEM velocity is equal to the 
required value, because the last point on the cr(V) profile, (Ttaem > is calculated from the 
QEGC by using the required values of r and V at TAEM. 

The variation of the velocity at the TAEM is monotonic with respect to o mid because 
larger value of cr mid means larger o(V) at any given V, hence smaller velocity at TAEM (due 
to faster energy dissipation); the opposite is also true. This monotonic functional behavior 
renders that the search for (7 mid to achieve required TAEM velocity will converge reliably 
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and quickly with a simple secant algorithm: Let Vtaem be the required TAEM velocity, and 
v taem the velocity at TAEM in the i-th search, corresponding to Then the secant 

algorithm is 


r *+ 1 

mid 


&mid 


&mid 


— a: 


i - 1 
mid 


Vtaem Vtaem 


{Vtaem ~ Vtaem) 


( 12 ) 


Note for each cr' mid a bank angle profile a{V) is defined, and Eq. (11) is integrated to get 
Vtaem • this process, the path constraints are enforced by imposing on a{V) the condition 

Vmin < <T(V) < (TmaxiV) (13) 


where <Tmax{V) is found in Eq. (9), and a min > 0 is a constant if an equilibrium glide 
condition at that constant bank angle is also desired to be enforced as in some cases . 4 The 
r- and V-versus-range profiles obtained will ensure that all the path constraints, range-to- 
HAC condition, and TAEM conditions on r and V are met. The QEGC not only enables 
us to reduce the need to numerically integrate more equations thus reduce the computation 
time, more importantly it allows us to enforce the path constraints so reliably and effortlessly. 


Once the above one-parameter search is completed, the obtained r and V profiles are 
appended to these obtained for the initial descent. Now we have r and V versus range-to-go 
profiles from entry to TAEM, and the magnitude of a along these profiles. 


2.4 Completion of 3D OF Trajectory 

Thus far the longitudinal state r and V versus range-to-go, and the magnitude of the 
bank angle have been determined. To complete the 3DOF trajectory design, the sign of 
a and the rest of the state need to be determined. This task is accomplished by tracking 
the obtained r and V versus range-to-go profiles with linear time-varying control laws for 
a and a ? From the r and V profiles, the corresponding 7 is easily backed out from Eq. 
(1). Note that the longitudinal dynamics Eqs. (1), (4) and (5) are decoupled from the rest 
of the state. Use the just obtained r, V, 7 , a (magnitude), and nominal a profile as the 
references. x; on — (r V 'y) T and u = (a a) T . Linearize the longitudinal dynamics, and still 
use range-to-go Stoga as the independent variable. Let the linearized longitudinal dynamics 
be 


Xlon A(sj 0 g 0 )(5xj 0n -f- B^Stog 0 )6\l 


(14) 


The trajectory control law developed in Ref . 2 will lead to 


— IC {3tago)^^-lon (^ 3 ) 

where the feedback gain K{st 0 go) ' s a function of A and B in Eq. (14) in closed-form. The 
modified bank angle a and angle of attack will be 


U = U Te f + <5u 


(16) 
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where u ref represent the bank angle profile obtained in preceding sections and nominal a. 
These controls in (16) are used to numerically integrate the complete 3DOF nonlinear dy- 
namics (1-6) from the entry interface. The range-to-go used to schedule the control law (15) 
is the actual range-to-go found from the integrated values of longitude 0 and latitude <j>. The 
initial sign of a is set to be the opposite of the sign of tx^HAC = (V> ~ ^wac)i where ^ hac 
is the azimuth angle to the HAC. The sign of o is reversed at a range s rev that is to be 
determined iteratively. The criterion for determining s rev is that A iJjhac at TAEM interface 
falls within the acceptable range (e.g., ±5 deg). Again the secant method appears to work 
well for this search. This search involves integration of all the dynamic equations (1-6), but 
does not need to be from the very beginning (entry interface). 

During this process, any limits on the magnitude of variations, rate and acceleration of 
the controls <r and ol can be applied. Once the integration (and search for s rev ) is completed, 
the 3DOF state trajectory and corresponding bank angle and angle of attack are obtained 
and stored as the reference trajectory. 

Compared to current standards in off-line trajectory planning, the above described tra- 
jectory generation algorithms produce a feasible trajectory very rapidly (about 2 seconds 
from a cold start, on a PC with a 800 MHZ processor). There are no apparent benefits in 
entry flight to require the generation of the reference trajectory in every guidance cycle. In- 
stead, The generated reference trajectory will be tracked by another guidance law similar to 
(15), but is based on the linearized dynamics of the overall system, not just the longitudinal 
dynamics. This guidance law will track both longitudinal and lateral state variables, and 
generate the actual guidance commands as described in Ref. 2 3 If necessary, the trajectory 
generation algorithms may be called again when the deviations from the reference trajectory 
become too large. In such a case the current condition of the RLV is used as the initial con- 
dition, and a new feasible reference trajectory is generated starting from the current state. 


4. Summary of the Algorithms 

The following is a re-cap of the steps described in the preceding sections for the trajectory 
generation algorithms: 

1. An initial descent phase when entry from orbit till sufficient dynamic pressure has been 
built up for transition to the QEGC. A constant bank angle is determined in this phase 
to ensure desirable descent behavior. The transition point is also found which marks 
the end of the initial descent phase. 

2. A single par am eter search defining the bank angle profile is done by using the QEGC 
and numerically integrating the velocity equation. This process determines the altitude- 
and velocity-versus-range profiles that satisfy all the path constraints and TAEM con- 
ditions on range-to-HAC, altitude and velocity. The magnitude of the bank angle along 
the QEGC is also found. 

3. Reduced order linear time-varying control laws for bank angle and ex are employed 
to track the just obtained altitude and velocity profiles. The complete set of 3DOF 
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dynamic equations are numerically integrated using the control law outputs. The 
sign of the bank angle is determined by a single bank reversal. The range at which 
the bank reversal occurs is iteratively determined by a secant method to reduce the 
heading-to-HAC error at the TAEM interface within the targeting limit. The final 
3D OF state trajectory and corresponding bank angle and angle of attack profiles are 
stored on-board as the reference trajectory. 

4. Linear time-varying trajectory control laws for o and ct based on the full-order 3DOF 
RLV dynamics are used on-line to track the reference trajectory. The actual guidance 
commands are the combinations of the reference bank angle and angle of attack and 
the outputs of these control laws. 

5. When necessary, new reference trajectory may be generated mid-stream by calling the 
algorithms again. 
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